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AN EMPIRICAL TECHNIQUE FOR COMPUTING NAVIGATION 
L COVARIANCE MATRICE8 AND THE APPLICATION OF THE TECHNIQUE TO 

POST APOLLO 8 LUNAR ORBIT MISSION PLANNING 
By John Bruce Williamson 


SUMMARY 




A n empirical technique used to calculate navigation covariance 
mu-ti ices is presented which is based on actual postflight tracking 
d ita processing from the Apollo 8 mission and the Langley Lunar 
i-oiter III. Postflight analyses of the Apollo 8 mission and 
Lunar Orbit er III revealed that the lunar gravitational potential 
function was not known as well as had been assumed before these missions. 
There is not sufficient information available at this writing to 
provide a satisfactory statistical description of the uncertainty in 
the model. This technique requires as input empirical errors in 
certain orbital parameters determined from postflight analysis. A 
numerical description (covariance matrix) of bhe anticipated navigation 
accuracies in lunar orbit is produced by the use of this technique. 

The technique is reasonably independent of the potential model which 
is used to process tracking data, although different values for the 
input parameters must be determined for each model. 

This technique has been applied to Apollo Missions F and G 
and has been found to model secular growth in down-track position 
errors satisfactorily. The technique is limited in its capability to 
lcdei the smaller growth rates expected in radial position errors, 
and it does not model oscillatory errors in the radial and down-track 
components of position and velocity as well as might be desired. A 
conservative estimate of the navigation accuracy is obtained when 
the oscillatory behavior is modeled as an independent secular error 
the magnitude of which is the amplitude of the oscillation. By careful 
application, t\ 4 s technique can be used to convert the results of 
postflight studies into the numerical input required for other mission 
planning tools (e.g., dispersion analysis programs). 
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INTRODUCTION 


Analysis of tracking from the Langley Lunar Orbiter program showed 
an unusual Doppler data residual pattern. Therefore, navigation of a 
spacecraft in lunar orbit with a perilune altitude of less then a 
few hundred miles could be inaccurate. Expex cs generally agree that 
this navigation problem is the result of an inaccurate lunar gravitational 
potential function. In other words, a significant amount of uneven 
mass distribution is present near the lunar surface. Analytic proce- 
dures for estimation of navigation accuracies require that all sources 
of error be precisely defined and that the probability distribution 
functions of the parameters be known. Studies are now being conducted 
to find the best lunar gravitational potential model for spacecraft 
navigation and the appropriate values to be assigned to the coefficients 
in the model. The probability distributions for the model are unknown 
at present. When a satisfactory model is implemented, the parameters 
of its probability distribution will be determined by processing 
tracking data from a number of different spacecraft. The primary 
effect of this problem is on processing earth-based radar dfta. 
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MSFN 
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PTEAP 

LOI-1 

LOI-2 
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RSS 
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ABBREVIATIONS AND SYMBOLS 

semimajor axis of spacecraft orb't 
Manned Spacecraft Center 
Manned Space Flight Network 

number of revolutions over which the covariance 
matrix is to be valid 

Perturbed Trajectory Error Analysis Program 

first lunar orbit insertion maneuver 

circularization maneuver 

spacecraft position vector 

root sum square 

Real-Time Computer Complex 

radius to spacecraft, r = |r| 
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s speed of spacecraft, s * |v| 

u v w coordinate spacecraft local orbit plane coordinate system defined 
system inert! ally at a g.ven time by the unit vectors 

i * R/r, J* (RxV)xR/| (R xV)xR| , k = RxV/|RxV| 


V 



Au , Av , etc . 


spacecraft velocity vector 


growth rates per revolution of the errors in local 
radius, down-track radius, orbital plane wedge anlge, 
and flight-path angle 


initial errors at the MSFN tracking epoch in 
local radius, down-track position, orbit plane 
wedge angle, and flight-path angle 


incremental change in u, v, etc. 


u 


gravitational constant of central body 


P 


U 



the correlation coefficient between the 
i^ h and parameters 


standard deviations in the components of the 
spacecraft position and velocity vectors, 
expressed in the u v w coordinate system 


a- 

w 


J 


u 


a , 

uv 


o . 
vu 


o 

Y 


covariance between radial position and 
down-track velocity 

covariance between down-track position 
and radial velocity 

standard deviation in flight-path angle 


o 


T 


standard deviation in orbital period 


T 


period of spacecraft orbit 


PROCEDURE 


The Mathematical Physics Branch of MSC has suspended the use 
of ordinary analytic tools for analyzing lunar orbit navigation 
based on the MSFN S-band radar stations. Instead, an interim 
technique has been devised which empirically models the trajectory 
effects observed in postflight processing of Lunar Orbiter III and 
Apollo 8 tracking data. This technique deals with trajectory parameters 
only and does not attempt to define the lunar potential model. Post- 
flight data processing provides the inputs to the technique as local 
errors and propagation effects observed with a particular lunar 
potential model. The resultant navigation accuracies are checked by 
propagating them with the standard analytic tools and by comparing 
these results with the postflight results. This assures that the 
application of the results (e.g. , dispersion analysis) using existing 
analytic tools will approximate the actual behavor of the spacecraft 
as observed in real time, to within the accuracy provided by the 
calibration with postflight results. 


ANALYSIS 


The current technique used to generate covariance matrices for 
vehicles in lunar orbit is based on the assumption that the navigation 
errors consist of independent local fit errors in four parameters and, 
independently, growth rates in these errors because of propagation. 
These parameters are as follows. 

u the local radius 

v the local down-track position 

0 the angle between the estimated and actual angular momentum 
vectors 

> the local flight-path angle of the vehicle 
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All other errors in the state vector of the vehicle are derived 
from this set. In general, the growth rate 0^ in 8 errors is so 

small that it is neglected for these studies. 

The Mr'FN tracking epoch is assumed to be in the middle of a 
front-side pass at approximately 0° longitude. This choice is 
arbitrary, and was made because it resulted in a formulation for 
out-of-plane errors which is convenient for lunar landing and rendezvous 
studies. If the tracking epoch is chosen to be in the middle of the 
data arc for one pass, the following equations are assumed for cross- 
track position and velocity errors. 

o * 0 

W (1) 

o. *= s sin (9 ) 
w u 

To account for some uncertainty in the longitude of the ascending 
node axid argument of latitude at the fit epoch is customarily 

set equal to 1000 feet. This does not result in a significant change 


The errors in radial position and flight-path angle must be 
constructed to include the error growth rate for a certain number 
of revolutions n. This construction is required because ordinary 
covariance propagation techniques do not model growths in these errors, 
which is one of the main restrictions on the technique. 


a 2 = (u n ) 2 + (n * u F 
u 0 g 

V = (y o )2 + (n * Y g )J 


(2) 


The growth rate in down-track position error is the most significant 
error growth, and this growth can be modeled with current propagation 
techniques. Therefore, the error in down-track position at the tracking 
epoch is simply 

o = v (3) 

v 0 


The development of the equations for the radial and down-track 
velocity errors and their covariances with the radial and down-track 
position errors are included in the appendix. In summary, the equations 
are the following. 
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(fc) 



All other correlation coefficients are zero. 

The covariance matrix at the MSFN tracking epoch may be found by 
use of equation (l) (modified) and equations (2) through (7). Because 
this matrix is symmetric, the numbers in the lower triangular portion 
ere replaced by their corresponding correlation coefficients, which 
nvey more intuitive information than covariances. The equation for 
t.i • correlation coefficient is 

p iJ * Jill (8) 


The covariance matrix is always presented in the vehicle's local 
u v w coordinate system defined inertially at the time associated 
with the matrix. 


RESULTS 


The error models empirically determined from unpublished post- 
flight tracking data analyses of the Langley Lunar Orbiter III and 
Apollo 8 are presented in table I. The three cases were used to provide 
navigation accuracies for Mission F planning and dispersion analysis 
support. Case 1 represents the assumptions that one pass of MSFN data 
will be processed in the RTCC orbit determination program with a 
triaxial lunar gravitational potential model and that both RTCC and on- 
board state vector predictors will use the triaxial model. Case 2 
represents the assumptions that one pass of MSFN data will be processed 
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with the Boeing R-2 lunar gravitational potential model (ref. l) 
and that both ground and onboard state vector predictors will use the 
R-2 potential model. Case 3 is the same as case 2 except that data 
from two consecutive passes will be processed. The Perturbed 
Trajectory Error Analysis Program (PTEAP) (ref. 2) was used to predict 
the st ;e vector uncertainties for the next two revolutions. The 
initial covariance matrices for these three cases are presented in 
table II. Plots of the predicted component and RSS position and velocity 
uncertainties for cases 1 through 3 are presented in figures 1 , 2, and 3. 

For comparison, postflight state vector comparisons from Apollo 8 
that correspond to case 3 are presented in figure U. A state vector 
was determined from the two-pass data arc in revolutions three and four 
by constraining the orbital plane to coincide with the orbital plane 
determined by a prediction of the pre-LOI-1 state vector through the 
confirmed LOI-1 and LOI-2 maneuvers. This vector was compared with 
the local one-pass unconstrained state vector solutions for revolutions 
four, five, and six. These comparisons are presented as differences in 
the components of the state vectors in the local u v w coordinate 
system. The view periods for this mission last from ho minutes after 
the beginning of the revolution until 1 hour 50 minutes after the 
beginning of the revolution. Comparisons between error analysis and 
postflight results should be made only during the view period. 


CONCLUSIONS 


The results presented are not intended to be independent deter- 
minations of navigation accuracies for Apollo missions. Rather, the 
technique is presented by which insight from pcstflight tracking data 
analyses is transformed into numerical input for the various studies 
which are performed in planning an Apollo-type mission. The figures 
should be compared with the postflight analyses to verify the technique. 
This comparison has been made informally and indicates that the 
technique models cases 1 and 2 quite accurately. In case 3, postf light 
analyses have revealed undamped (i.e., growing in amplitude) oscillatory 
errors, especially in radial position and down-track velocity (fig. U), 
which are not modeled by this technique. However, the technique does 
model the secular effects satisfactorily. A conservative estimation 
of the navigation accuracies is obtained if the oscillatory behavior 
is modeled as an independent secular error the magnitude of which 
is the amplitude of the oscillation. Case 3 is by far the most 
satisfactory from a navigation standpoint. The cross-track errors 
in case 3 are the same as in case 2 because the cross-track components 
of the state vector are determined by a one-pass data fit, and the two- 
pass fit is constrained to maintain these values. Recent unpublished 
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nostflight analyses Indicate that a better estimation of cross-track 
position and velocity during lunar orbit operations can be obtained by 
prediction of the pre-!iOI-l state vector through the confirmed maneuvers. 
This technique would be used with the one-pass technique as a monitoring 
tool. Hence, the error estimates described in this note should be 
conservative in cross-track position and velocity. 




.Parameter 

Case I 

Case II 

Case III 

u , ft 

0 

1000 

1000 

1000 

v , ft 

6 

3000 

3000 

3000 

• 

• 

• 

• 

• 

• 

• 

• 

<u 

Ti 

0k 

O 

CD 

0.1 

0.1 

0.1 

Y . deg 

.01 

.01 

.01 

u ft /rev 

g 

1500 

500 

500 

v , ft/rev 

g 

15 000 

11 000 

1500 

0 , deg/rev 

.005 

.005 

.005 

Y * deg/rev 

6 

.005 

.005 

.005 
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TABLE II.- COVARIANCE MATRICES FOR LUNAR ORBIT 


MANNED SPACE FLIGHT NETWORK (MSFN) TRACKING DATA 


(a) Navigation covariance matrix for case 1, one-pass triaxial fit 

at MSFN tracking epoch a 


9.9999996 * .10 6 0 0 

0 9.0 x 10 6 0 

0 0 1.0 x 

0 -.78U66998 0 

-.97036671 0 Q 


0 -8809.^81+3 0 

-7928.5361 0 0 

10 6 0 00 
11.31+1+081+ 0 0 

0 8.21+2951+7 0 


0 


0 0 0 


0 



(b) Navigation covariance matrix for case 2, 


one-pass R-2 fit at 


2.0 x 10 C 0 0 

0 9.0 x 10 6 0 

0 0 1.0 x 10 6 

0 -.891+51+89 0 

-.92565272 0 0 

0 0 0 


MSFN tracking 

epoch a 


0 

-1761.8969 

0 

-7928.5361 

0 

0 

0 

0 

0 

8.7281+11+1 

0 

0 

0 

I.811I+85I+ 

0 

0 

0 

88.063182 


£ 

Format for matrices in table II. 
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TABLE II.- COVARIANCE MATRICES FOR LUNAR ORBIT 

MANNED SPACE FLIGHT NETWORK (MSFN) TRACKING DATA - Continued 

(c) Navigation covariance matrix for case 3, 
two-pass R-2 fit, at MSFN tracking epoch 8. 


2.0 x 10 
0 
0 
0 

-0. 9981*5 
0 


9.0 x io 
0 

- 0 . 89U55 
0 
0 


0 

0 

1.0 X 10* 
0 
0 
0 


-7928.536 

0 

8.7281*12 

0 

0 


Format for matrices in table II. 
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(a ) 2 
u 

0 

0 

0 
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uv 
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uv 


(a ) 2 
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vu 
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0 


0 

( °v ): 
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0 

0 


vu 

0 


V 2 


0 

0 


(a.) 2 


-1761.897 

0 

0 


1.556963 

0 
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0 

0 

0 

0 


<°w > 2 


0 

0 

0 

0 

0 

81.0 
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TABLE II.- COVARIANCE MATRICES FOR LUNAR ORBIT 
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(a) One sigma position. 

Figure 1.- One pass MSFN solution accuracies, triaxial lunar 

potential model. 
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Time elapsed from MSFN tracking epoch, min 


(b) One sigma velocity. 
Figure 1 . - Concluded. 






(a) One sigma position. 

Figure 2.- One pass MSFN solution accuracies, Boeing R-2 lunar 

potential model. 


1 tr velocity uncertainty, fps 


l6 
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(b) One sigma velocity. 
Figure 2 . - Concluded. 





Figure 3.- Two pass MSFN solution accuracies, Boeing R-2 lunar potential model with 
out-of-plane solution constrained to one pass R-2 solution. 
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(b) One sigma veloc.ty. 







16:00 


Greenwich mean time, G.m.t., hrrmin 


(a) Radial position difference. 


Figure 4.- Apollo 8 state vector from revolutions 3 and 4, constrained fit compared 
with local one revolution unconstrained fits. 




































Revolution 4 


Revolution 5 


Revolution 


Greenwich mean time, G.m.t., hrrmin 


(d) Radial velocity difference 








Revolution 6 


Revolution 4 


Revolution 5 


Greenwich mean time, G.m.t., hr:min 


(e) Downtrack velocity difference 


Figure 4.- Continued 
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(f) Crosstrack velocity difference. 


Figure 4 . - Concluded . 
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APPENDIX 

DERIVATION OF THE EQUATIONS FOR THE ERRORS 
IN RADIAL AND DOWN-TRACK VELOCITY AND THEIR COVARIANCES 
WITH RADIAL AND DOWN-TRACK POSITION ERRORS 
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APPENDIX 


DERIVATION OF THE EQUATIONS FOR THE ERRORS 
IN RADIAL AND DOWN-TRACK VELOCITY AND THEIR COVARIANCES 
WITH RADIAL AND DOWN-TRACK POSITION ERRORS 

RADIAL VELOCITY 


The radial velocity error is assumed to be a function of flight 
path angle error and phasing (i.e., down-track position) errors. 
Assuming that these two error sources are independent, one may write 



(Al) 



The equation relating radial velocity to flight-path angle is 


u = s sin y 


so that 
orbits , 



we may 


cos y . 
write 


Since we are dealing with near-circular 


Y - 0 

cos y - 1 (A2) 

9u * s 

3y 

For circular orbits, in which the flight-path angle is constant, a 
phasing error translates into an error in the direction of the vfiocity 
vector, of the same magnitude. Thus an angular phasing error <t> 
translates into a radial velocity error 



\ 

, * 

; / 
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Au = -s sin 4 > 


where 


So 


or 


sin 4> 




Av 

r 


Au = — Av 
r 


3u s_ 

3v r ( A3 ) 


Substituting equations (A2) and (A3) into (Al) we have 



(AU) 


DOWN-TRACK VELOCITY 


The down-track velocity error is also assumed to be caused by 
independent error sources. The first of these is the relationship 
between speed and radius derived from the energy equation assuming 
a fixed semimajor axis. The second is a period error. The period 
error translates directly into an error in semimajor axis and the 
energy equation is used to relate this error to speed. This results 
in the following equation 



Taking the partial derivatives 
elliptical orbits 


of the energy equation for 



s 


z 


V 
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We have 


27a 2 


*. 


(A6) 


© 


The relationship between period and semima.lor axis is obtained from 


t = 2 n — 

A U 


which can be solved for semimajor axis 


O 


l/2 T 


which differentiates to yield 


Ba _ 2a 
3t * 3t 


( AT ) 


For nearly circular orbits the following relat ionsnips hold 

u = r and v - s 


and from the energy equation 


a - r and u - rv^ . 


(A8) 
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An additional assumption involved in this derivation is that the 
down-track position error growth rate is a result of the period error 
which is contained entirely in the down-track velocity error at the 
t.‘*"Oking epoch. This assumption can be stated as 



(A9) 


Substituting the relationships (A8) into equations (A6) and (AT), and 
substituting this result along i/ith ( A9 ) into equation (A5) we have 



(A10) 


COVARIANCES 


Since the uncertainties in the in-plane velocity components 
were computed from the uncertainties in the in-plane position 
components, these errors are correlated. We may compute the correlation 
coefficients as follows 


c . 
uv 




(All) 


c 


vu 



£ 

r 


o 

v 


2 


All other correlation coefficients are zero. 
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